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ABSTRACT 



The interstellar abundances of D I, N I, and O I in the local ISM are stud- 
ied using high-resolution spectra of four hot white dwarfs. The spectra of GD 
246, WD 2331-475, HZ 21, and Lan 23 were obtained with the Far Ultravio- 
let Spectroscopic Explorer (FUSE) in the wavelength range 905 - 1187 A. The 
line of sight to GD 246 probes the Local Interstellar Cloud and at least one 
other H I cloud inside the Local Bubble, which contains most of the gas seen 
along this line of sight. The column densities of H I, C II*, S II, and Si II are 
measured using archival Hubble Space Telescope STIS echelle-mode observations. 
The H I column density is determined by fitting the strong damping wings of 
interstellar Lya using a model atmosphere to account for the stellar continuum. 
The sightline-averaged ratios for GD 246 are: D I/H I = (1.51 ± x 10-^ 
O I/H I = (3.63 ± I'll) X 10-^ and D I/O I = (4.17 ± x 10"^ (uncertainties 
are 2a). This line of sight provides the fourth reliable FUSE measurement of the 
Local Bubble D/H ratio. For the WD 2331—475 line of sight we find sightline- 
averaged ratios: D I/O I = (5.13 ± jfg) x 10"^ and D I/N I = (4.57 ± HI) x 10"^ 
Toward HZ 21 the sightline-averaged ratios are: D I/O I = (4.57 ± I'H) x 10"^ 
and D I/N I = (4.27 ± x 10^^. In the higher column density sightline to 
Lan 23 the sightline-averaged ratios are: D I/O I = (3.24 ± ^Qg) x 10"^ and 
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D I/N I = (3.16 ± X 10~^. Molecular hydrogen, corresponding to rotational 
levels J < 3, is clearly seen along this line of sight. No reliable H I measurements 
are available for WD 2331—475, HZ 21, or Lan 23. We combine the different 
abundance ratios computed here with previous published values to produce re- 
vised FUSE abundance ratios for D I/H I, O I/H I, N I/H I, D I/N I, D I/O I, 
and O I/N I. 

Subject headings: ISM: Abundances — ISM: Evolution — Ultraviolet: ISM — 
Stars: Individual (GD 246, WD 2331-475, HZ 21, Lan 23) 

1. INTRODUCTION 

The present day abundance ratio of deuterium to hydrogen places important constraints 
on Big Bang nucleosynthesis (BBN) and the chemical evolution of galaxies. Since it is be- 
lieved that deuterium is only produced in appreciable amounts in primordial BBN and 
destroyed in stellar interiors, the measurement of D I/H I in the interstellar medium (ISM) 
places a lower limit on the primordial abundance of deuterium. By comparing the ISM 
abundance of deuterium to its abundance in high-redshift intergalactic gas we should be 
able to understand better the effects of astration and chemical evolution of galaxies. Mea- 
surements of the D/H ratio in intervening clouds of gas seen toward distant quasars have 
yielded a range of values D/H = (1.65 - 4.0) xlO"^ (O'Mcara et al. 2001; Pettini & Bowen 
2001; Levshakov et al. 2002, and references therein). Measurements of D/H in the local 
ISM have been made with Copernicus (e.g. Rogerson & York 1973), HST, (e.g. Linsky et al. 
1995), IMAPS, (Jenkins et al. 1999; Sonneborn et al. 2000), and more recently FUSE (Moos 
et al. 2002, and references therein). A nearly constant ratio of D/H = (1.5 ± 0.1) x 10^^ 
{la on the mean) has been obtained in the Local Interstellar Cloud (LlC) by Linsky (1998); 
recent measurements inside the Local Bubble (Moos et al. 2002) appear to be consistent with 
a single value for D/H in the Local Bubble. Other measurements (Laurent et al. 1979; York 
1983; Hebrard et al. 1999; Jenkins et al. 1999; Sonneborn et al. 2000) suggest variations of 
the interstellar D/H ratio beyond the Local Bubble, at the distance of a few hundred parsecs. 

Until FUSE was launched (1999 June 24) only the UV spectrographs onboard HST could 
be used for systematic studies of the deuterium abundance in the ISM. Since only the 
Lya transitions of H I and D I are observable in the HST bandpass, such studies were 
restricted to low column density sightlines so that D I was neither completely obscured by 
the adjacent H I Lya nor heavily saturated. With FUSE, we now have access to the com- 
plete Lyman series of deuterium (except Lya), which allows the study of a larger range of 
environments and provides tighter constraints on the abundance of D I compared to studies 
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using only Lya. In addition, the numerous lines of N I, O I, and Fe II in the FUSE bandpass 
can be used to trace the metallicity and dust content of the absorbing gas. 

In this work we present the first measurements of deuterium absorption toward the 
white dwarfs GD 246, WD 2331-475, HZ 21, and Lan 23. Table 1 summarizes the stellar 
parameters. This study adds four more sightlines to the seven sighthne studies that constitute 
the first set of FUSE deuterium measurements (see overview paper by Moos ct al. 2002, and 
references therein). This paper is organized as follows. The observations and data processing 
are presented in §2; the analysis methodology is described in §3. The GD 246 line of sight 
analysis is presented in §4. The WD 2331—475, HZ 21, and Lan 23 lines of sight are analyzed 
in §5, and systematic effects are discussed in §6. We discuss the results in §7. 



2. OBSERVATIONS AND DATA PROCESSING 

2.1. FUSE Observations 

The FUSE observatory consists of four co-ahgned prime-focus telescopes and Rowland- 
circle spectrographs that produce spectra over the wavelength range 905 - 1187 A with a 
spectral resolution of ~ 15 - 20 km s~^ (wavelength dependent), for point sources. Two 
of the optical channels employ SiG coatings, providing reflectivity in the wavelength range 
~ 905 - 1000 A, while the other two have LiF coatings for maximum sensitivity above 1000 A. 
Dispersed light is focused onto two photon-counting microchannel plate detectors. With this 
arrangement of optical channels (LiF 1, LiF 2, SiC 1, and SiG 2) and detector segments (lA, 
IB, 2A, 2B) the FUSE instrument has 8 segments: LiF lA, LiF IB, LiF 2A, LiF 2B, SiG lA, 
SiG IB, SiG 2A, and SiG 2B. Four channels cover the wavelength range 1000 - 1080 A while 
two channels each cover the ranges 900 - 1000 A and 1080 - 1180 A. The FUSE mission, its 
planning, and on-orbit performance are discussed by Moos et al. (2000) and Sahnow et al. 
(2000) 

Table 2 summarizes the FUSE observations of the four white dwarfs studied in this work. 
The FUSE spectrum of the hot DA white dwarf GD 246 is shown in Figure 1. Observations 
were obtained through both the large (LWRS, 30" x 30") and the medium sized apertures 
(MDRS, 4" X 20") in histogram mode (HIST, in which a two-dimensional spectral image 
is accumulated and downlinked). The spectrum of WD 2331—475 is shown in Figure 2. 
Two observations were obtained through the LWRS aperture and one through the MDRS 
aperture. All of them were obtained in HIST mode. During one of the LWRS observations 
(P1044202) one detector (side 2) was not available resulting in the loss of data from the LiF 
2A, LiF 2B, SiG 2A, and SiG 2B channels. The spectrum of HZ 21 is shown in Figure 3. The 
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data were obtained through the LWRS and MDRS apertures, in time-tag mode (TTAG, in 
which the position and arrival time of each photon is recorded). The FUSE spectrum of the 
hot DA white dwarf Lan 23 is displayed in Figure 4. The data were obtained through the 
LWRS aperture in TTAG mode. 

For the four stars the two-dimensional FUSE spectra are reduced using the CalFUSE 
pipeline (version 2.0.5, Dixon, Kruk, & Murphy, in preparation^). The processing includes 
data screening for low quality or unreliable data, thermal drift correction, geometric distor- 
tion correction, heliocentric velocity correction, dead time correction, wavelength calibration, 
detection and removal of event bursts, background subtraction, and astigmatism correction. 
The spectra are aligned by cross-correlating the individual exposures over a short wavelength 
range that contains prominent spectral features and then coadded by weighting each expo- 
sure by its exposure time. All the spectra are binned to three pixel samples, or 20 mA, for 
analysis (the line spread function, LSF, is about 11 pixels or ~ 70 mA wide). 



2.2. STIS Observations of GD 246 and Data Processing 

Two HST /STIS observations of GD 246 that include H I Lya arc available at the 
Multimission Archive at the Space Telescope Science Institute; one each using the E140H 
and E140M gratings. 

Although the E140H grating has better resolution than the E140M grating it also suffers 
from more complications in the Lya region. Lya spans more orders in the E140H grating 
and the relative flux cahbration of adjacent orders is not as good as the one for the E140M 
grating. Because the saturated H I Lya hue is so broad, the higher resolution of the E140H 
grating is not required. Therefore we make use of the E140M data for deriving the H I 
column density. Wc use the E140H data to study the abundance of metals such as C II*, Si 
II, and S II and to determine the velocity structure of the gas along the line of sight. 

The STIS echelle-mode observations of GD 246 are summarized in Table 3. The E140M 
(S/N ~ 15, FWHM = 1.3 pix) and E140H (S/N ~ 10 - 25, FWHM = 1.2 pix) data are 
reduced using the standard STSDAS pipeline within IRAF. The CALSTIS pipeline (ver- 
sion 2.3) is used to produce the two-dimensional spectrum of the FUV-MAMA and then 
the scattered-light removal is performed using different algorithms for the different echelle 
gratings. For the E140M data the Lindler & Bowers (2000) algorithm is used to estimate 
and remove the scattered light from the data in the pipeline procedure. After extracting the 



^The CalFUSE pipeline reference guide is available at littp://fuse. plia.jliu.edu/analysis/pipeline_reference. html 
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spectrum, two echelle spectral orders in the regions adjacent to Lya are combined using a 
weighted averaging scheme where the orders overlap. For the E140H data the scattered light 
is removed with the procedure of Howk & Sembach (2000). 

3. Analysis Methodology 

Whenever possible we derive the column densities by more than one method. We use 
profile fitting (PF), curve of growth (COG) and apparent optical depth (AOD) methods. 
Table 4 indicates which methods arc used, for each transition, in the four stars. Below we 
discuss the techniques used in more detail. 

3.1. Profile Fitting 

We use the profile fitting code Owens. f developed by Martin Lemoine and the French 
FUSE Team to derive the column densities presented in this study. In this approach each 
interstellar absorption line is represented by the convolution of a theoretical Voigt absorp- 
tion profile with the instrument line spread function (LSF), taken to be a single Gaussian 
in this work. An iterative procedure which minimizes the sum of the squared differences 
between model profiles and the data is used to determine the most likely column densities, 
N, radial velocities, and velocity dispersions of each component. The stellar continua in the 
vicinity of the absorption lines being studied are normalized by Ghebyshev polynomials of 
low order. Unsaturated lines from different elements in different spectral regions (designated 
"windows") are fit simultaneously to give a single solution for the temperature, non-thermal 
velocity, and radial velocity of each cloud (we define as unsaturated, those lines for which the 
absorption profile, before convolution with the instrumental LSF, has a residual intensity in 
the line core larger than 0.1). The software also allows us to use a wavelength dependent LSF 
and relative wavelength shifts between different absorption lines as free parameters. These 
last two features are particularly convenient when fitting FUSE data, for which the LSF, 
while poorly constrained, is known to vary with wavelength. There are also small wavelength 
shifts between and within channels which are accommodated by this approach. For more on 
profile fitting with this code see Hebrard et al. (2002) and Lemoine et al. (2002) 

Two of us (G. Oliveira and G. Hebrard) performed the analyses for GD 246, WD 
2331—475, HZ 21, and Lan 23 independently. In both of the analyses a Gaussian LSF 
is used. However, in one analysis (PFl) the FWHM of the LSF for each spectral window 
is allowed to vary independently during the fitting process, while in the other (PF2) it is 
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not allowed to vary. In PF2 there is an initial step in which all LSF's are free parameters 
of the fit. After a few trial fits the LSF's for the different windows converge and do not 
change significantly (these values are close to each other and close to what is known about 
the FUSE LSF). At this point the LSFs are fixed at the newly found values, reducing the 
number of free parameters, for the remainder of the minimization. We find typical LSF 
values of 8 - 13 pixels (FHWM, 53.6 - 87.1 mA) for the PFl analyses; for the PF2 analyses 
we find a range of 9 - 13 pixels for the LSF (FWHM, 60.3 - 87.1 mA). The background levels 
used for each fine are those evaluated at the bottom of the closest Lyman fine (see Hebrard 
et al. 2002). Where possible, we compare the different channels and datasets to check for 
fixed-pattern noise in the regions of the absorption lines used in this study. When there is 
an inconsistency between channels and/or observations we remove that absorption line from 
our fit. 

The error bars for each estimated column density are computed from an analysis of the 
variation. The column density for which the error bars are being sought is fixed at a series 
of trial values, and for each trial value we compute the best fit and resulting distribution 
where all the other parameters are allowed to vary freely. To define the confidence interval 
we use the distribution of x^- Scanning the column density in this way, we obtain the value 
of x^ as a function of the column density and derive the 1, 2, 3, 4, 5, 6, and 7a confidence 
levels using the standard A^^ method. The 1, 2, 3, 4, 5, 6, and 7a error bars are divided 
by 1, 2, 3, 4, 5, 6, and 7, respectively. The average of these values, multiphed by two, is 
then used to produce the final 2a uncertainties reported here. For more on the Ax^ method 
see Hebrard et al. (2002) or Lemoine et al. (2002). We adopt the wavelengths and oscillator 
strengths for molecular hydrogen from Abgrall et al. (1993a,b) and use Morton (1991) with 
updates from Morton (1999, private communication, see Table 4) for the other species. 

Both analyses produce similar results and the final values quoted here reflect the com- 
bined effort of these analyses. The results for all the four lines of sight are summarized in 
Table 5 below, in § 5. 

3.2. Curve of Growth and Apparent Optical Depth 

We use the curve of growth (COG) and apparent optical depth (AOD) techniques to 
estimate the column densities of some species as a consistency check of the proflle fltting 
results. Not all species analyzed by proflle fltting are suitable for analysis with these two 
methods. In the cases where the continuum is difficult to place (such as D I), or where 
there is blending (such as the N 1 triplets) we do not construct a curve of growth, or use 
the apparent optical depth technique. Since in these two techniques each absorption line 
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is studied individually, continuum placement in a "problematic" line is harder then when 
many different lines of different absorption strengths and different continua are analyzed 
simultaneously as is the case of profile fitting. For species with only a small number of 
transitions available (such as P 11 and Ar I) we do not construct a curve of growth but we 
are able to use the apparent optical depth method. 

In the COG technique the measured equivalent widths are fit with a single component 
Gaussian curve of growth (Spitzer 1978). In this technique the Doppler parameter, b, and 
column density, A^, are varied so as to minimize the between the measured equivalent 
widths and a model curve of growth. The stellar continuum in the vicinity of each line is 
estimated using a low order Lcgcndrc polynomial fit to the data. Contributions from Poisson 
noise, uncorrected high frequency fixed pattern noise, uncertainties in the Legendre fit pa- 
rameters, and systematic uncertainties in the continuum placement and velocity integration 
range are used to estimate the uncertainties in the measured equivalent widths (see Sembach 
& Savage 1992, for details on this procedure). We measure the equivalent widths of a par- 
ticular transition in all segments where that wavelength is covered. All the measurements 
from the different segments are compared; those that differ by more than 2a are excluded 
from our analysis (the differences are probably due to fixed pattern noise). 

In the AOD technique the column density is determined by directly integrating the 
apparent column density profile, Na{v) = 3.768 x 10^'^Ta{v)/[fX{A)], over the velocity 
range of the absorption profile (sec Savage & Sembach 1991). The continuum is placed in 
the same way as described above in the COG method. This technique only yields the true 
total column density when the absorption is weak (r < 1) or if partially saturated, the 
components of the fines are fully resolved. 

All the lines used in PF, COG, and AOD analyses arc shown in Table 4 along with 
log /A for each transition. For each star, P, C, and A denote transitions that are used in 
profile fitting, curve of growth, and apparent optical depth, respectively. Both the curve of 
growth and apparent optical depth techniques produce results which generally have higher 
uncertainties, but are accurate enough to use as consistency checks. These checks show that 
the AOD and COG results are within the 2a uncertainties of the column densities determined 
using profile fitting. 

4. ANALYSIS OF THE GD 246 LINE OF SIGHT 

GD 246 is a hot hydrogen-rich (DA) white dwarf with effective temperature and gravity 
Te// ~ 53,000K and log g ~ 7.85. The FUSE spectra of GD 246 shows a smooth stellar 
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continuum against which the interstellar lines are cleanly seen. Apart from the broad stellar 
Lyman lines, only a few stellar lines are present (Si IV, P IV-V). No molecular hydrogen 
and O VI lines are detected along this line of sight. 

High-resolution STIS echelle spectra of GD 246 show that there are at least two velocity 
components, separated by ~ 9 km s~^, along this hne of sight. These are seen in the N I, S 
II, and Si II absorption profiles (see Figure 5). As discussed in § 7.5, the weaker component 
likely probes the LIC. 

The two profile fitting analyses of this line of sight are different in an additional way. In 

one of the profile fitting analyses (PFl) of this line of sight we use only FUSE data with a 
variable LSF; in the other profile fitting analyses (PF2) wc use both FUSE and STIS data 
with fixed, but different LSFs, for the FUSE and STIS spectral windows. Both analyses are 
described in more detail below. 



4.1. Profile Fits 

In PFl wc fit a one component model to D I, C II*, N I, O I, Si II, P II, Ar I, and 
Fe II using FUSE data only. The FHWM of the LSF is a free parameter of the fit for each 
spectral window being fit. 

For PF2 we use high resolution STIS E140H spectra (R ~ 114,000) of GD 246, covering 
the wavelength region 1170 A - 1372 A, to constrain the velocity structure along the line 
of sight as well as to measure the column density of species with no transitions in the 
FUSE bandpass (such as S II). Species which have transitions in the FUSE and STIS 
bandpasscs arc fit simultaneously in both datascts. For the STIS data wc use a single- 
Gaussian LSF with a FWHM of 2 pixels, for all the spectral windows. In this analysis we 
fit one absorption component to D I, C II*, N I, O I, Si II, S II, and Fe II, at the best-fit 
velocity of Vq — —9.76 km s~^. For N I, Si II, and S II we use two components; component 
1 at the best-fit velocity of Vq — —9.76 km s~^ and component 2 at the best-fit velocity of 
Vq — —0.63 km s"-*^. 

The decision to distribute the second set of ions among 2 components is based on the 
following arguments. From an independent analysis of the STIS data alone, where we fit 
two components to N I, Si II, and S II, we determined that component 1 contains most 
of the material seen along this line of sight and that the Doppler parameters of the two 
components are similar. The ratio of column densities of component 1 to component 2 are 
N(N I)i/N (N 1)2 = 5.9, N(S II)i/N(S 11)2 = 4.6, and N(Si II)i/N (Si 11)2 = 2.3. For 
species for which there is no STIS data, a velocity separation of ~ 9 km s"^ between the two 
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components is not enough for profiles at the FUSE resolution to produce meaningful results 
about the relative strength of the two different components. Since the physical parameters 
of the components seem to be similar we would expect component 1, the stronger one, to 
saturate first. However, we chose lines in the FUSE data that appear unsaturated so that 
we did not have to worry about hidden saturation effects that would bias our column density 
determinations. For these reasons we choose to fit two components for the species for which 
we have STIS echelle data and one component for the species with no STIS data. In the 
end, the results of the two different approaches, PFl and PF2, are consistent at the 2a level, 
with most consistent at the Icr level. 

Because of the better resolution and signal-to-noise ratio of the STIS E140H data, we 
have tighter constraints on the column densities of C II*, S II and. Si II, in the model where 
STIS and FUSE data are fitted simultaneously. For this reason the column densities and 
error bars reported for these species (Table 5) are the ones obtained with PF2. Profile fits 
for D I (PF2) are presented in Figure 6 and the combined results are presented in Table 5. 

4.2. Determination of the H I Column Density Toward GD 246 

The H I column density along this sightline is determined by fitting the strong damping 
wings of interstellar Lya using a newly computed model atmosphere to account for the 
stellar continuum. We consider the infiuence on the Lya profile of several stellar models 
with different effective temperatures and gravities. ifST/STIS E140M (A/AA = R 45,800) 
archival data are used in this analysis. 

4.2.1. Stellar Model 

Several investigators measured the atmospheric parameters of GD 246 by fitting the 
Balmer lines with the use of either LTE or non-LTE (NLTE) stellar atmospheres models. 
They obtained effective temperatures and gravities that range from 53,000 K < T^q < 60,000 
K and 7.7 < log^ < 8.0 (see, e.g., Finley et al. 1997; Marsh et al. 1997; Venues et al. 1997; 
Barstow et al. 1998; Napiwotzki et al. 1999). For our analysis of the stellar continuum in the 
vicinity of the H I Lya profile, and knowing that NLTE effects are important in the core of 
Lya, we adopted the atmospheric parameters determined by Napiwotzki et al. (1999) who 
based their results on pure-hydrogen NLTE models. They obtained T^s — 53,088 ± 968 K 
and loggr = 7.85 ±0.07 (Icr error bars). Napiwotzki et al. (1999) pointed out that systematic 
uncertainties, which may result from the extraction of the spectrum, the fiux calibration, or 
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the normalization, are much larger than the statistical errors. To take into account these 
systematic errors, we considered the la errors on T^s and logg to be 3,000 K and 0.15. 
The grid of models that we computed to estimate the stellar Lja line includes these errors 
(Teff = 50,000 K, 53,000 K, and 56,000 K; log 5 = 7.70, 7.85, and 8.0). 

Observations in the EUV, FUV, and UV wavelength bands reveal that the atmosphere of 
GD 246 contains low abundances of elements heavier than hydrogen. For instance, Chayer 
et al. (2001) and Wolff et al. (2001) analyzed the International Ultraviolet Explorer and 
FUSE spectra of GD 246, and measured very low abundances of carbon, silicon, and phos- 
phorus. Wolff et al. (1998) analyzed the Extreme Ultraviolet Explorer spectrum of GD 246 
with LTE metal-line blanketed models and concluded that the star had a lower mctallic- 
ity than G191-B2B, which is the metal-rich white dwarf prototype. Dupuis et al. (2000) 
reported the detection of Fe VI absorption features in the CHANDRA LETG spectrum of 
GD 246. But the non-detection of Fe V in the STIS spectrum of GD 246 puts an upper limit 
on the iron abundance of log(Fe/H) < —6.0. We expect that these small traces of heavy ele- 
ments do not modify the physical structure of the atmosphere of the star, and consequently, 
do not influence the Lja line profile. Therefore, we computed pure-hydrogen NLTE stellar 
atmospheres models to describe the stellar Lya line. 

Determining the H I column density with different stellar models places more credible 
error bars on N(H I) than using a single stellar model. The statistical uncertainties are 
derived by using the best stellar model, and the systematic uncertainties associated to the 
stellar models arc estimated by considering the most extreme stellar parameters, i.e., those 
parameters that yield the strongest and weakest stellar Lja. Figure 7 illustrates such stellar 
Lya profiles. The figure shows that the best stellar model (model a; T^s = 53,000 K and 
logs' — 7.85) produces a Lya line profile that is between the strongest Lya line (model 
b; Teff = 50,000 K and log^ = 7.70) and the weakest one (model c; T^s = 56,000 K and 
logy = 8.00). 

4.2.2. Lya Profile Fitting 

Wolff et al. (2001) determined that the photospheric lines of this star are separated by 
V — —6 km s~^ from the interstellar fines. To take into account the relative shift between 
interstellar and stellar lines we align the model and the data in three different ways; displacing 
the model by —10, 0, and -|-10 km relative to the data. By comparing the stellar Lya 
model to the data, we took into consideration the radial velocity of GD 246, which includes its 
velocity with respect to the Sun and its gravitational redshift. Because of the large breadth 
of the interstellar H I absorption, displacing the model and the data by 10 km s~^ does not 
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lead to any difference in tfie measured H I column density from our fits. For each stellar 
model computed, the E140M spectrum of H I Lja is normalized by the stellar model prior 
to profile fitting (see Figure 7, bottom panel). In addition, a 2"^^ order polynomial is used 
to model the continuum after normalization by the model atmosphere, to take into account 
uncertainties in the instrument sensitivity. Using the best fit stellar model we perform two 
fits (see § 3 for a description of the fitting process) of the H I Lya profile using both a 2 
pixel single-Gaussian LSF and the tabulated STIS LSF for the 0'.'02 x (y.'09 sht. There are no 
differences in the results using the two LSFs. We derive log N(H I) = 19.11 ± 0.02 ± 0.03 
(2(7) where the first error bars reflect the uncertainties associated with the different stellar 
models and the second error bars take into account the statistical uncertainties when using 
the best fit model. Wc find then, log N(H 1) = 19.11 ± 0.05 {2a), by taking a conservative 
approach to combine the two different uncertainties. Our H 1 column density is in close 
agreement with the value of log N(H1) = 19.12 ± q^q (lo") derived from EUVE data by 
Barstow et al. (1997) using stellar models incorporating heavy elements. 

Lemoine et al. (2002) have expressed concern over the presence of weak H 1 components 
which could perturb or bias the measurement of N(H 1), but which would not be detected 
in other species due to their weak column density. Vidal-Madjar & Ferlet (2002) reanalyzed 
the Capella line of sight and showed a significant increase in the uncertainty due to this 
effect. However, they also concluded that for warm gas with column densities of H I above 
10^^ cm~^ the H I column density evaluation is more credible than for absorbers with lower 
column densities. We have tried to investigate the possible impact of additional weak, hot 
components on the determination of N(H I). In order to do so, we added one or two additional 
components in H I only and performed the fits to Lya and to Ly/?. For Ly/? we used the 
blue wing of night only data, since the geocoronal emission in this region distorts the shape 
of the red wing of the profile. We were not able to find weak hot components of the type 
described by Lemoine et al. (2002) for which both Lya and Ly/5 profiles agreed. As pointed 
out by Lemoine et al. (2002) these weak hot components are very difficult to find and it is 
possible that we might have missed them, thus underestimating the uncertainties quoted for 
N(H I). However, we think this is unhkely given the conclusions of Vidal-Madjar & Ferlet 
(2002). 

5. ANALYSIS OF THE WD 2331-475, HZ 21, AND Lan 23 LINES OF 

SIGHT 

WD 2331-475 is a hot DA white dwarf with %e = 51,800 K and log g = 7.79 (Vennes 
et al. 1997). A few stellar lines of Si IV, P VI-V, S IV-V are present in the data. There are 
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no HST observations of this star and the quahty of the lUE data is not high enough for a 
rehable H I measurement. 

HZ 21 is a hehum rich (DO) white dwarf showing only a few stellar hues of N IV, P V, Si 
IV and S VI in the FUSE data. Figure 3 shows the broad stellar He II stellar absorption hues 
due to transitions from the energy level n — 2. With effective temperature Teff = 53,000 K 
and gravity log g — 7.8 (Dreizler & Werner 1996) this low metallicity star provides a smooth 
stellar continuum against which the interstellar lines are clearly seen. No molecular hydrogen 
is detected along this line of sight. There is probably more than one velocity component 
along this line of sight, considering the large distance to this star, = 115 pc. Unfortunately 
no HST observations of this star with sufficient resolution to resolve the velocity structure 
exist. In addition, this star has not been observed by EUVE and the quality of the lUE data 
is not high enough for a reliable H I measurement. 

Lan 23 is a DA white dwarf with an effective temperature of T^s = 59,360 K and log 
g — 7.84 (Venues et al. 1997). Stellar lines from the ground states of O VI, S VI, and 
from an excited state of C HI are seen in the FUSE spectra of this star. This star has 
not been observed by HST, so we do not have information about the velocity structure 
along this hue of sight. Considering that this star is at a distance oi d — 122 pc, several 
absorption components would be expected. Wolff et al. (1999) have used EUVE data to 
measure the interstellar H I column density along this line of sight. Because of the high 
interstellar column density of H I and also some unknown opacity in the white dwarf itself, 
the EUVE fiux is very small, which implies a large uncertainty in the determination of 
the interstellar H I column density. They find log N(H I) = 19.89 ± 004, where the error 
bars quoted include uncertainties in the photospheric composition. These authors do not 
mention the confidence level corresponding to the quoted error bars, so we assume these are 
approximately la uncertainties, which we use to calculate the 2a uncertainties quoted in 
Table 7. 

5.1. Column Density Determinations 

5.1.1. WD 2331-475 

For this line of sight we measure the column densities of D I, C II*, N I, O I, Si II, P II, Ar 
I, and Fe II, using profile fitting, COG and AOD. Both PFl and PF2 discussed previously 
were used. We first attempted to fit the data with a single absorption component. This 
resulted in a large underestimation of the D I, N I, and O I column densities when compared 
with the COG and AOD results. A closer inspection of the data revealed that some of the 
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stronger lines have a wing extending towards negative velocities. Profile fitting with two 
components resulted in column density values within la of the COG and AOD results. We 
determined that the data is better represented by two absorption components separated by 
~14 km s~^, with similar Doppler parameters, b. Figure 8 displays the two component fit 
to C II*, N I, and Fe II. This component structure is also used for D I, O I, Si II, P II, and 
Ar I, not shown in Figure 8. The fit to D I (using PF2) is shown in Figure 9. The column 
densities adopted are given in Table 5. 

5.1.2. HZ 21 andLan 23 

For both of these sightlincs our profile fitting analyses assume the presence of only one 
interstellar component. We fit the species D I, N I, O I, Si II, P II, Ar I, and Fe II in HZ 
21 and D I, N I, O I, and Fe II in Lan 23. For the Lan 23 line of sight we include an extra 
component containing H2, described in more detail in §5.2. 

In the case of HZ 21, the uncertainties introduced by the presence of broad stellar He II 
absorption in the neighborhood of the D I lines should be negligible compared to our final 
error bars on the column density of D I. Dreizler & Werner (1996) determined an upper limit 
on the hydrogen to helium abundance of H/ He < 0.1 in the atmosphere of the helium-rich 
white dwarf HZ 21. To take into consideration the effects of the stellar hehum abundance 
on the column density of D I, we computed two stellar models with hydrogen to helium 
abundances of H/He = 0.1 and 0.01. Figure 10 displays two stellar models computed with 
Tcff = 53,000 K and log (7 = 7.8, but different H to He abundance ratios. Model a has 
H/He = 0.1 and model h has H/He = 0.01. In both models the He II stellar absorption is so 
broad that it does not affect the underlying D I absorption. 

There are many oxygen transitions in the FUSE bandpass covering a wide range of 
oscillator strengths. Unfortunately at the high column density of O I in Lan 23 and with 
the low signal-to-noise ratio of the data, it is difficult to find non-saturated O I lines that 
are suitable for our analysis. We find only two such lines, A919.917 and A974.070, which 
span more than 1 dex in log/ A. This is the source of the large error bars for the O I column 
density along this line of sight. 

Fits for D I (using PF2) are presented in Figures 11 (HZ 21) and 12 (Lan 23). The 
column densities found by combining the results of the two different analyses (PFl and PF2), 
are given in Table 5. 
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5.2. Interstellcir H2 TowEird Lan 23 

Absorption arising from H2 rotational levels J = - 3 is clearly seen in the FUSE spec- 
tra of Lan 23. Some H2 lines are blended with D I and can be a problem for the measurement 
of the D I column density. We determine the column density of each J level by fitting a 
single component curve of growth to lines from J = to 3, simultaneously, assuming the 
same h value for all J levels. Prom the curve of growth we obtain a Doppler parameter, 
6 = 4.1 ± km s~^ (2(j). The h obtained from the curve of growth is then used as a 
constraint in the profile fitting analysis of H2, described in the paragraph below. 

We also determine the H2 column density by including unsaturated H2 J = - 3 lines 
in our profile fitting analysis of this line of sight. The column density for each rotational 
level is determined independently, with no assumption about the rotational temperature of 
the gas. To allow for the possibility that this small amount of molecular hydrogen does not 
reside in the same gas as the D I, the velocity centroid of the H2 is free to assume a different 
value from that of the D I. 

Column densities obtained with curve of growth and profile fitting are displayed in 
Table 6. The total H2 column density along this line of sight is of the order log N(H2) ~ 15. 
The derived excitation temperatures are Tqi ~ 291 K, T12 ~ 324 K, and T23 ~ 240 K, 
however, formal relative contributions of UV pumping versus coUisional excitation have not 
been explored. 

6. SYSTEMATIC EFFECTS 

There are several systematic effects in the FUSE data which can contribute to uncertain- 
ties in the estimated column densities and that are not directly included in our determination 
of the error bars. Below we discuss the most important systematic effects in detail, how they 
can affect our data, and how we address these problems in our profile fitting analysis. For 
additional discussion of systematic effects in D/H measurements such as this, we refer the 
reader to Hebrard et al. (2002); Kruk et al. (2002); Moos et al. (2002). 

6.1. Fixed-Pattern Noise 

Fixed-pattern noise associated with the detectors can distort profile shapes and alter 
the equivalent widths of the absorption lines. The four channel design of FUSE provides at 
least two (sometimes four) independent spectra of most of the lines used in the analysis. By 
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comparing data from multiple segments, multiple observations and multiple slits we are able 
to identify the absorption profiles that suffer from large (compared to the statistical noise) 
fixed-pattern noise artifacts and exclude them from our fits. However, small effects cannot 
be ruled out. 



6.2. Background and Line Spread Function 

Although a comprehensive study of the line spread functions for the FUSE instrument 
is not yet available, preliminary studies indicate that it is composed of two (Gaussian) 
components. The narrow component has a FWHM of ~ 15 - 20 km s~^ and the broad 
component contains 20 - 30% of the total LSF area (with a large portion of this broad 
component falling under the narrow core of the LSF). This broad component is responsible 
for the residual light in the cores of saturated lines (e.g., the high-order Lyman lines in 
Figures 1-4); in addition there may be scattered hght. We use a single Gaussian LSF in 
our analysis and adopt a small zero flux level shift to account for the residual light in lines 
close to the Lyct series of H I, where this effect is more visible. Because we are fitting only 
non-saturated lines, any residual errors in this correction for the tail of the LSF should 
not affect our derived column densities in a significant way. To the extent that this small 
amount of residual light is due to the LSF, it may produce small systematic changes in the 
derived column densities, but such effects are expected to be small compared to the derived 
uncertainties. 



6.3. Oscillator Strength Uncertainties 

While oscillator strengths (/-values) of the D I lines are well known, the same may not 
be true for other elements such as N I and O 1. By fitting multiple lines of the same species 
at the same time we are able to remove the lines that stand out as a poor fit and average 
out the small errors of the other lines in the final result. A case where the uncertainty in an 
O I /-value might pose a problem is in the determination of the O I column density for the 
Lan 23 line of sight. For this star, only two non-saturated lines have enough signal-to-noise 
ratio to be useful in our analysis: A974.07 (very weak, / = 1.56x10"^ and partially blended 
with H2(J = 2)) and A919.9 (/ = 1.76x10-^). 

In other analyses of O I where A919.9 was used (e.g. Hebrard et al. 2002) this fine did 
not stand out as having a large problem with its /-value. A test, where the O I column 
density was determined by fitting each line individually, showed a 20% dispersion around 
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the value obtained with a fit to all the fines simultaneously (see Hebrard et al. 2002). Since 
the error bars derived for N(0 I) toward Lan 23 are large, it is unlikely that the uncertainty 
in the /-values of the O I lines used here play a large role in our measurement of N(0 I). 

We measure the column density of Ar I in the four lines of sight using only one transition 
of Ar I (A1048.2198) because the other one (A1066.6599) is blended with photospheric Si IV. 
Any uncertainties in the /-value of this transition will affect our column density of Ar I. 

6.4. Number of Interstellar Components Along the Line of Sight 

In § 4.2.2 we discussed the potential effects of weak hot H I components on the deter- 
mination of N(H I) for GD 246. Because the H I column density is greater than 10^^ cm~^, 
such effects are probably small. For all species a major difficulty associated with measuring 
column densities is the presence of iinrcsolvcd overlapping absorption components along the 
line of sight. Small scale velocity structure has been seen in, amongst others, K 1 and Na 
I ultra-high resolution surveys, for relatively small path lengths (see for example Welty & 
Hobbs 2001). The problem arises if the physical conditions in the different absorption com- 
ponents are significantly different. In a mix of cold (narrow) and warm (broad) components, 
each with approximately the same column density, the cold component will become satu- 
rated before the warm component. Because the components are unresolved, the absorption 
profile resulting from the sum of the two components may appear to be unsaturated, thus 
leading to an incorrect estimate of the total column density. For GD 246 the STIS data 
shows no evidence of such material. Lan 23 has a low H2 column density (~ 10^^ cm~^), 
and the other white dwarfs show no evidence of H2 in their FUSE spectra. However, without 
ultra-high resolution data this systematic effect can never be completely ruled out. 

6.5. Profile Contamination by Unknown Stellcir Features 

Although improbable, some of the lines used in our analysis could be contaminated by 
unknown stellar features (e.g. Sonncborn et al. 2002; Friedman et al. 2002). Since the four 
white dwarfs studied here are all metal poor, we do not believe this is a significant concern. 
Photospheric Si IV is present in the spectra of the white dwarfs studied here and is blended 
with Ar I A1066.66; this Ar I transition is not used in our analysis. 
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6.6. Effect on Column Density Ratios 

We believe that the systematic effects discussed above are generally small. It is possible 
that these effects could be of importance for the column densities with small statistical 
uncertainties. However, it is unlikely that they will be of significance for the column density 
ratios discussed in 5 7. 



7. RESULTS AND DISCUSSION 

In section 7.1 we discuss the line of sight to the four stars in the context of the Local 
Bubble. In § 7.2 we compare the total column density ratios of D/H, 0/H, and N/H for 
GD 246, WD 2331-475, HZ 21, and Lan 23 to values found for other Unes of sight in the 
local ISM and discuss the imphcations of using O and N as tracers of H for these lines of 
sight. In § 7.3 we combine the ratios derived here with previously published FUSE ratios, 
deriving revised FUSE mean ratios. In § 7.4 we discuss the mean D/0 ratio inside the 
Local Bubble. § 7.5 presents a more detailed discussion of the GD 246 line of sight. Because 
we have information on the velocity structure along this line of sight we can separate the 
Lie and Local Bubble contributions to the total column densities and draw more detailed 
conclusions about the distribution of material along the line of sight. 

7.1. Local Bubble 

Sfeir et al. (1999) have used absorption line studies of Na I to approximately map 
the contours of the Local Bubble. They found a cavity with a radius between 65 and 250 
pc (depending on the direction) that is delineated by a sharp gradient in the neutral gas 
column density with increasing radius, a dense neutral gas "wall". Sfeir et al. (1999) quote 
an equivalent width of 20 mA of Na I A5891.59 as corresponding to log N(H I) ~ 19.3, which 
is used as a rough estimate of a plausible Local Bubble boundary. According to their map 
of the Local Bubble, the GD 246 hne of sight does not penetrate the wall, as it lies between 
the 5 mA and 20 mA Na I A5889.95 contours. Although the uncertainty in the distance to 
GD 246 (and for the other three white dwarfs, for that matter) makes this result ambiguous, 
log N(H I) = 19.11 ± 0.05 implying that GD 246 is hkely close to the wall of the Local 
Bubble, but does not penetrate the higher density regions beyond. The lines of sight to 
WD 2331-475 and HZ 21 do not penetrate the wall either; WD 2331-475 lies inside the 20 
mA contour and HZ 21 lies between the 5 mA and 10 mA contours. In addition, as shown in 
the next section, the measured O I column density toward both WD 2331—475 and HZ 21 
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implies an H I column density similar to that towards GD 246. The Lan 23 line of sight is 
located between the 20 mA and 50 mA Na I A5889.95 contours, and has a high H I column 
density; it likely penetrates the wall of the Local Bubble. 

7.2. O/H, D/H, and N/H Ratios 

Table 7 summarizes the column density ratios for several species along the four lines 
of sight studied in this work. For GD 246 we find O 1/ H I = (3.63 ± g:^^)xlO-^. This 
ratio is in agreement with the O I/H I ratio derived by Meyer et al. (1998), of O I/H 
I = (3.43 ± 0.30) X 10"^ (using the /-value updated to / = 1.16x10"^ from Welty et al. 
(1999). 

Rehable H I measurements are not available for WD 2331—475, HZ 21, or Lan 23. 
However, the constancy of the O/H ratio discussed above enables us to derive an approximate 
column density for H 1 using O I. Using the Meyer et al. (1998) interstellar O 1 abundance 
and the O I column density along each line of sight, we estimate log N(H I) = 18.94 for WD 
2331-475, log N(H 1) = 19.20 for HZ 21, and log N(H I) = 20.16 for Lan 23. 

Using the N(H I) obtained this way and our measured N(D I) we obtain D 1/H 1 ~ 1.8x10"^, 
1.6xl0-^ and 1.2xl0^^ for WD 2331-475, HZ 21, and Lan 23, respectively The uncer- 
tainties in the O 1 columns used to derive the H I column densities are large enough that 
the differences between the derived D/H ratios for these stars and the Moos et al. (2002) 
weighted mean of D/H = (1.52 ± 0.15) xlO"^ are not significant. The D/H ratio for GD 246 
is discussed below in § 7.5. 

Using the Meyer et al. (1997) interstellar abundance of nitrogen, N/H = (7.5 ± 0.8) xlO~^, 
and our measured N I column density leads to estimated values of log N(H I) = 18.87 for 
GD 246, log N(H I) = 18.65 for WD 2331-475, log N(H I) = 18.89 for HZ 21, and log N(H 
I) = 19.82 for Lan 23. The H I column density derived using the N I column density is in all 
four cases only ~50% of the H I column density determined using the O I abundance, imply- 
ing that the N I column densities are low. It is possible that the N along these lines of sight is 
~50% ionized. This effect has been seen for other sightlines as well (see Jenkins et al. 2000; 
Lehner et al. 2002; Moos et al. 2002). For GD 246 we measure N 1/H 1 = (4.37 ± ^of^)xlO-^ 
which is a factor of 1.7 times smaller than the Meyer et al. (1997) value, but is in close 
agreement with the weighted mean of the FUSE sighthnes N I/H I = (4.24 ± 0.62) xlO~^ 
(see Moos et al. 2002). 

The ionization fraction of N I along a line of sight can be estimated directly by mea- 
suring the N II and N HI column densities as both N II and N HI have transitions in the 
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FUSE bandpass. However, for our sightlines the N II transitions available in the FUSE band- 
pass are saturated, and N III cannot be determined since it is blended with interstellar Si 
II and sometimes with photospheric N III. Hebrard et al. (2002); Kruk et al. (2002); Wood 
et al. (2002) have been able to measure the column density of N II (and sometimes N III) 
along three different lines of sight and found that in all cases it was larger than the column 
density of N I, providing a better agreement between {N/}i)totai within 100 pc and N 1/B.totai 
for larger path lengths. 

7.3. Revised FUSE Ratios 

In Table 8 we combine the ratios derived here with the Moos et al. (2002) ratios for 
D I/H I, O I/H I, N I/H I, D I/N I, D I/O I, and O I/N I, and present the Meyer et al. 
(1997, 1998) N/H and 0/H values for comparison. We present also the xl test for a single 
mean value. The values quoted in Table 8 are the weighted means and uncertainty in the 
means, where we use the largest of the lower and upper error bars of each individual ratio 
to compute the weighted mean. We use la error bars in this subsection in order to make 
it easier to compare the different lines of sight among themselves. For Lan 23 values only 
the D I/N I value is used, due to the large error bars for the other ratios associated with 
this line of sight . By combining our ratios with the Moos et al. (2002) values, the weighted 
means increase by to 5%. 

The mean D I/O I ratio is (4.06 ± 0.17) xlO'^ with xl = 1-9 (z/ = 9). The mean D I/N 
I ratio now computed with ten lines of sight is (3.41 ± 0.15) xlO~^ with xl — 2.2 {u — 9). 

7.4. D/O Ratio Inside the Local Bubble 

Hebrard et al. (2001) and Moos et al. (2002) have found a nearly constant D/O ratio 
inside the Local Bubble, indicating that D/O traces the D/H ratio. The weighted mean of the 
five FUSE sightlines within the Local Bubble (Moos et al. 2002) is (3.76 ± 0.20) xlO^^ 
the mean). When we combine the D/O ratios derived here (Lan 23 value not included, as it is 
likely outside the Local Bubble) with the Local Bubble FUSE results from Moos et al. (2002), 
we derive a weighted mean and uncertainty in the mean of D I/O I = (3.87 ± 0.18) x 10^^ (la 
in the mean) with xl = 0.96 = 7), reinforcing the Moos et al. (2002) conclusion that the 
variability of the data is consistent with the uncertainties, i. e., the D I/O I ratio is constant 
within the Local Bubble. The D I/O I ratio inside the Local Bubble will be discussed in 
more detail by Hebrard et al. (2002, in preparation). 
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7.5. GD 246 

Absorption by at least two different components, separated by ~9 km s~^, is clearly 
seen, along the GD 246 sightline, in the high resolution STIS echelle data of N I, S II, and Si 
II (see Figure 5). The stronger component (component 1), seen at —9.76 km s^^, contains 
86%, 82%, and 79% of the total column densities of N I, S II, and Si II, respectively, along 
this hne of sight. The velocity of component 2, V2 — —0.63 km s~^, is consistent with the 
expected velocity of the Local Interstellar Cloud (LIC) along this hne of sight (+2.46 ± 1.84 
km s^^; Lallement et al. 1995) when the zero point velocity uncertainty in STIS E140H data 
is taken into account {Av ~ 0.7 - 1.5 km s~^). The discrepancy between the expected LIC 
velocity and the measured velocity along this sightline might be an indication that the LIC 
velocity field is more complicated than a single vector. A similar effect has been seen toward 
the Hyades Cluster by Redfield & Linsky (2001), who found f predicted (LIC) — t'observed(LIC) 
= 2.9 ± 0.7 km s~^, similar to our value Vpredicted(LIC) — Vobserved(LIC) ~ 3.1 km s~^. The 
eighteen Hyades stars studied by these authors are located in the sky not too far frm GD 
246; for four of them (HD 27561, HD 27848, HD 28237, and HD 29225; I ~ 178 - 184, b ~ 
-20 25) they also find Vedicted(LIC) - t'observed(LIC) ~ 3.0 km s~^. 

Sulfur should not be depleted onto dust grains, making it possible to estimate how 
neutral hydrogen would be distributed between the two components using S II as a proxy 
for H I. We derive log N(H 1)2 = 18.39 using the S II column density for the LIC component 
(component 2), log N(S 11)2 — 13.59. We derive a distance of 8 pc to the edge of the LIC 
along this direction (assuming a LIC neutral hydrogen column density of Uhi = 0.10 cm~^; 
Linsky et al. 2000). Both the H I column density and distance to the edge of the LIC 
derived here are greater than the values derived from the model of Redfield (2002, private 
communication), for this direction, of log N(H I) = 17.93 and d = 2.75 pc. Our values are 
closer to the maximum values derived for the LIC by Redfield & Linsky (2000) of 18.32 and 
6.8 pc, respectively. The discrepancy between our values and the LIC model values can be 
due either to the preliminary nature of their model or to the presence of an unresolved cloud 
close to the LIC velocity. While H I may not follow exactly the same distribution as S II 
it is still clear that most of the neutral gas along this line of sight is outside the LIC but 
within the Local Bubble boundaries. 

We find D I/H I = (1.51 ± g:i)xlO-^ for this line of sight {2a error bars). So far, 
only three lines of sight studied by FUSE (HZ 43A, G191-B2B, and WD 1634-573) probing 
the Local Bubble have accurate D/H measurements (Kruk et al. 2002; Lemoine et al. 2002; 
Wood et al. 2002). This is the fourth line of sight, with an accurate D/H measurement that 
probes the region mentioned above. 

The D I/H I ratio for aU the gas in this sightline is (1.51 ± gJ^xlO-^ Linsky (1998) 
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reports a mean value of (1.50 ± 0.10) xlO~^ (Icr in the mean) for 12 LIC sightlines. Because 
most of the gas along the sightline to GD 246 resides outside the LIC, we conclude that the 
D I/H I ratio for the gas outside, is also ~1.5xl0~^, similar to the LIC value. 

The electron density in component 1, which probes material outside the LIC, can be 
computed under the assumption that electron collisions populate the upper fine structure 
level (J = 3/2) of the ground state of C II leading to C II* (see Spitzer & Fitzpatrick 1993). 
Because all the C II lines in the FUSE and STIS data are saturated, S II is used as a proxy 
for C II. This is a reasonable assumption; C II and S II have similar ionization potentials 
(24.38 eV for C II and 23.33 cV for S II) and should therefore trace the same type of gas. 
The coUisional equilibrium equation is 

A2in{C+*) = 7i2nen(C+) (1) 

where 712 is the rate coefficient for excitation and A21 is the spontaneous downward transition 
probability. For excitation of ions by electrons (Spitzer 1978): 

8.63x10-6^ , s -1 

712 = — ^^12 exp{-—) cm s (2) 

= 7.9x10-2 eV, and we adopt A21 = 2.29 (Mendoza 1983) and ^12 = 2.90 (Ostcrbrock 
1989). Assuming that conditions are uniform within the cloud we can replace particle density 
ratios by column density ratios and assuming that for warm low density clouds the depletions 
of S and C are and —0.3 dex respectively, wc find 

N{SII) ■ ^ ^ 

Using the S II column density we measured for component 1, log N(S II)i = 14.25 ± 0.02, and 
the C II* column density, log N(C II*) = 13.05 ± 0.04, we find Ug = (0.1 ± 0.01)(T/8000)0-5 
cm^"2 (all error bars are 2o"). For T = 4000-8000 K this yields an electron density, rie = 0.06- 
0.11 cm~^ for the material beyond the LIC, which is similar to values found in the LIC by 
Linsky et al. (2000). 

Thus, the D/H ratio and the electron density obtained in this study indicate that the 
properties of the cloud(s) beyond the LIC along the GD 246 line of sight are likely similar 
to those of the LIC. 

The presence of ionized gas in this cloud will not affect the D/H ratio in this region, as 
D I and H I, with a similar ionization potential, are ionized at the same rate. 

Table 9 compares abundances toward GD 246 with solar photosphere abundances. We 
find there is an indication that N, O, Si, P, Ar, and Fe are depleted from typical solar 
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abundances, but the uncertainties preclude a definitive statement. Moos et al. (2002) and 
Linsky (2002) have noted the evidence for 20 - 25% depletion of O I by dust in the near 
ISM. For all species listed above part of this depletion might be explained by not taking into 
account higher ionization states when calculating their gas-phase abundances. N and Si are 
the only species for which we could measure some of this contribution; however we have only 
a lower limit on N(N II) and Si III is blended with stellar Si III. 

7.6. Summary 

We have obtained column densities of H I, D I, C II*, N I, O I, Si II, P II, S II, Ar I, 
and Fe II for the line of sight to the white dwarf GD 246. This is the fourth line of sight 
with an accurate D/H measurement with FUSE that probes the region outside the LIC but 
inside the Local Bubble. Based on a constant D/H ratio of 1.5x10^^ inside the LIC we find 
a similar ratio outside the LIC. In addition, our sightline-averaged D I/H I measurement for 
GD 246 is in agreement with the weighted mean of D I/H I reported by Moos et al. (2002) 
of D I/H I = (1.52 ± 0.08)xl0-^ for the local ISM. For the WD 2331-475, HZ 21, and 
Lan 23 sightlines the column densities of D I, N I, O I, and other species are also measured 
and ratios are computed. We have used results from all four sightlines to compute revised 
FUSE ratios for the local ISM. 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by The Johns Hopkins University. Financial support to 
U. S. participants has been provided in part by NASA contract NAS5-32985 to Johns Hopkins 
University. Support for French participation in this study has been provided by CNES. Based 
on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data 
Archive at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with proposal 7296. The the profile fitting procedure, Owens.f, 
used in this work was developed by M. Lemoine and the French FUSE Team. 

REFERENCES 

Abgrall, H., Roueff, E., Launay, F., Roncin, J. Y., & Subtil, J. L. 1993a, A&AS, 101, 273 
— . 1993b, A&AS, 101, 323 

AUende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 556, L63 



-23- 



Barstow, M. A., Dobbie, P. D., Holberg, J. B., Hubeny, I., & Lanz, T. 1997, MNRAS, 286, 
58 

Barstow, M. A., Hubeny, I., & Holberg, J. B. 1998, MNRAS, 299, 520 

Chayer, R, Kruk, J. W., Vennes, S., & Dupuis, J. 2001, in ASP Conf. Ser. 226: 12th 
European Workshop on White Dwarfs, 90 — h 

Dreizler, S. & Werner, K. 1996, A&A, 314, 217 

Dupuis, J., Vennes, S., & Chayer, P. 2000, American Astronomical Society Meeting, 197, 
Finley, D. S., Koester, D., & Basri, G. 1997, ApJ, 488, 375 

Friedman, S. D., Howk, J. C, Chayer, P., Tripp, T. M., Hebrard, C, Andre, M., Oliveira, 

C, Jenkins, E. B., Moos, H. W., Oegerle, W. R., Sonneborn, C, Lamontagne, R., 
Sembach, K. R., & Vidal-Madjar, A. 2002, ApJS, 140, 37 

Grevesse, N. & Sauval, A. J. 1998, Space Science Reviews, 85, 161 

Hebrard, C, Friedman, S., Kruk, J. W., Lehner, N., Lcmoinc, M., Linsky, J. L., Moos, 
H. W., Oliveira, C. M., Sembach, K. R., Sonneborn, C, Vidal-Madjar, A., & Wood, 
B. E. 2001, in Gaseous Matter in Galaxies and Intergalactic Space, 18th lAP meeting 
held June 19-23, 2001, Paris. Ed. Ferlet et al., p. 79. 

Hebrard, G., Lemoine, M., Vidal-Madjar, A., Desert, J.-M., Lecavelier des Etangs, A., Ferlet, 
R., Wood, B. E., Linsky, J. L., Kruk, J. W., Chayer, P., Lacour, S., Blair, W. P., 
Friedman, S. D., Moos, H. W., Sembach, K. R., Sonneborn, G., Oegerle, W. R., & 
Jenkins, E. B. 2002, ApJS, 140, 103 

Hebrard, G., Mallouris, C, Ferlet, R., Koester, D., Lemoine, M., Vidal-Madjar, A., & York, 

D. 1999, A&A, 350, 643 

Holweger, H. 2001, in AIP Conf. Proc. 598: Joint SOHO/ACE workshop "Solar and Galactic 
Composition", 23 — h 

Howk, J. C. & Sembach, K. R. 2000, AJ, 119, 2481 

Jenkins, E. B., Oegerle, W. R., Cry, C, Vallcrga, J., Sembach, K. R., Shelton, R. L., Ferlet, 
R., Vidal-Madjar, A., York, D. G., Linsky, J. L., Roth, K. C, Dupree, A. K., & 
Edelstein, J. 2000, ApJ, 538, L81 



Jenkins, E. B., Tripp, T. M., Wozniak, P. . A., Sofia, U. J., & Sonneborn, G. 1999, ApJ, 
520, 182 



-24- 



Kruk, J. W., Howk, J. C, Andre, M., Moos, H. W., Oegerle, W. R., Oliveira, C, Sembach, 
K. R., Chayer, P., Linsky, J. L., Wood, B. E., Ferlet, R., Hebrard, G., Lemoine, M., 
Vidal-Madjar, A., & Sonneborn, G. 2002, ApJS, 140, 19 

Lallement, R., Ferlet, R., Lagrange, A. M., Lemoine, M., & Vidal-Madjar, A. 1995, A&A, 
304, 461 

Laurent, C, Vidal-Madjar, A., & York, D. G. 1979, ApJ, 229, 923 

Lehner, N., Gry, C, Sembach, K. R., Hebrard, G., Chayer, R, Moos, H. W., Howk, J. C., 
& Desert, J.-M. 2002, ApJS, 140, 81 

Lemoine, M., Vidal-Madjar, A., Hebrard, G., Desert, J.-M., Ferlet, R., Lecavelier des Etangs, 
A., Howk, J. C., Andre, M., Blair, W. P., Friedman, S. D., Kruk, J. W., Lacour, S., 
Moos, H. W., Sembach, K., Chayer, P., Jenkins, E. B., Koester, D., Linsky, J. L., 
Wood, B. E., Oegerle, W. R., Sonneborn, G., & York, D. G. 2002, ApJS, 140, 67 

Levshakov, S. A., Dessauges-Zavadsky, M., D'Odorico, S., & Molaro, P. 2002, ApJ, 565, 696 

Lindler, D. & Bowers, C. 2000, American Astronomical Society Meeting, 197, 

Linsky, J. L. 1998, Space Science Reviews, 84, 285 

— . 2002, Space Science Reviews, in press 

Linsky, J. L., Diplas, A., Wood, B. E., Brown, A., Ayres, T. R., & Savage, B. D. 1995, ApJ, 
451, 335 

Linsky, J. L., Redfield, S., Wood, B. E., & Piskunov, N. 2000, ApJ, 528, 756 

Marsh, M. C, Barstow, M. A., Buckley, D. A., Burleigh, M. R., Holberg, J. B., Koester, D., 
O'Donoghue, D., Penny, A. J., & Sansom, A. E. 1997, MNRAS, 286, 369 

Mendoza, C. 1983, in lAU Symp. 103: Planetary Nebulae, Vol. 103, 143-172 

Meyer, D. M., Cardelh, J. A., & Sofia, U. J. 1997, ApJ, 490, L103+ 

Meyer, D. M., Jura, M., & Cardelh, J. A. 1998, ApJ, 493, 222 

Moos, H. W., Cash, W. C, Cowie, L. L., Davidsen, A. F., Dupree, A. K., Feldman, P. D., 
Friedman, S. D., Green, J. C, Green, R. F., Gry, C, Hutchings, J. B., Jenkins, 
E. B., Linsky, J. L., Malina, R. F., Michalitsianos, A. G., Savage, B. D., ShuU, 
J. M., Siegmund, O. H. W., Snow, T. P., Sonneborn, G., Vidal-Madjar, A., Wilhs, 
A. J., Woodgate, B. E., York, D. G., Ake, T. B., Andersson, B.-G., Andrews, J. P., 



-25- 



Barkhouser, R. H., Bianchi, L., Blair, W. P., Brownsberger, K. R., Cha, A. N., 
Chayer, P., Conard, S. J., Fullcrton, A. W., Gaines, G. A., Grange, R., Gummin, 
M. A., Hebrard, G., Kriss, G. A., Kruk, J. W., Mark, D., McCarthy, D. K., Morbey, 

C. L., Murowinski, R., Murphy, E. M., Oegerle, W. R., Ohl, R. G., Ohveira, C., 
Osterman, S. N., Sahnow, D. J., Saisse, M., Sembach, K. R., Weaver, H. A., Welsh, 
B. Y., Wilkinson, E., & Zheng, W. 2000, ApJ, 538, LI 

Moos, H. W., Sembach, K. R., Vidal-Madjar, A., York, D. G., Friedman, S. D., Hebrard, 
G., Kruk, J. W., Lehner, N., Lemoine, M., Sonneborn, G., Wood, B. E., Ake, T. B., 
Andre, M., Blair, W. P., Chayer, P., Cry, C, Dupree, A. K., Ferlet, R., Feldman, 
P. D., Green, J. C, Howk, J. C, Hutchings, J. B., Jenkins, E. B., Linsky, J. L., 
Murphy, E. M., Oegerle, W. R., Oliveira, C, Roth, K., Sahnow, D. J., Savage, B. D., 
ShuU, J. M., Tripp, T. M., Weiler, E. J., Welsh, B. Y., Wilkinson, E., & Woodgate, 
B. E. 2002, ApJS, 140, 3 

Morton, D. C. 1991, ApJS, 77, 119 

Napiwotzki, R., Green, P. J., & Saffer, R. A. 1999, ApJ, 517, 399 

O'Meara, J. M., Tytler, D., Kirkman, D., Suzuki, N., Prochaska, J. X., Lubin, D., & Wolfe, 
A. M. 2001, ApJ, 552, 718 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galactic nuclei (Research 
supported by the University of California, John Simon Guggenheim Memorial Foun- 
dation, University of Minnesota, et al. Mill Valley, CA, University Science Books, 
1989, 422 p.) 

Pettini, M. & Bowen, D. V. 2001, ApJ, 560, 41 

Redfield, S. & Linsky, J. L. 2000, ApJ, 534, 825 

— . 2001, ApJ, 551, 413 

Rogerson, J. B. & York, D. G. 1973, ApJ, 186, L95+ 

Sahnow, D. J., Moos, H. W., Ake, T. B., Andersen, J., Andersson, B.-G., Andre, M., Artis, 

D. , Berman, A. F., Blair, W. P., Brownsberger, K. R., Calvani, H. M., Chayer, P., 
Conard, S. J., Feldman, P. D., Friedman, S. D., Fullerton, A. W., Gaines, G. A., 
Gawne, W. C, Green, J. C, Gummin, M. A., Jennings, T. B., Joyce, J. B., Kaiser, 
M. E., Kruk, J. W., Lindler, D. J., Massa, D., Murphy, E. M., Oegerle, W. R., Ohl, 
R. G., Roberts, B. A., Romelfanger, M. L., Roth, K. C, Sankrit, R., Sembach, K. R., 
Shelton, R. L., Siegmund, O. H. W., Silva, C. J., Sonneborn, G., Vaclavik, S. R., 
Weaver, H. A., & Wilkinson, E. 2000, ApJ, 538, L7 



-26- 



Savage, B. D. & Sembach, K. R. 1991, ApJ, 379, 245 

Sembach, K. R. & Savage, B. D. 1992, ApJS, 83, 147 

Sfeir, D. M., Lallement, R., Crifo, F., & Welsh, B. Y. 1999, A&A, 346, 785 

Sonneborn, G., Andre, M., Oliveira, C, Hebrard, G., Howk, J. C, Tripp, T. M., Chayer, P., 
Friedman, S. D., Kruk, J. W., Jenkins, E. B., Lemoine, M., Moos, H. W., Oegerle, 
W. R., Sembach, K. R., & Vidal-Madjar, A. 2002, ApJS, 140, 51 

Sonneborn, G., Tripp, T. M., Ferlet, R., Jenkins, E. B., Sofia, U. J., Vidal-Madjar, A., & 
Wozniak, R . R. 2000, ApJ, 545, 277 

Spitzer, L. 1978, Physical processes in the interstellar medium (New York Wiley-Interscience, 
1978. 333 p.) 

Spitzer, L. J. & Fitzpatrick, E. L. 1993, ApJ, 409, 299 

Vennes, S., ThejU, P. A., Galvan, R. G., & Dupuis, J. 1997, ApJ, 480, 714 

Vidal-Madjar, A. & Ferlet, R. 2002, ApJ, 571, L169 

Welty, D. E. & Hobbs, L. M. 2001, ApJS, 133, 345 

Welty, D. E., Hobbs, L. M., Lauroesch, J. T., Morton, D. C., Spitzer, L., & York, D. G. 
1999, ApJS, 124, 465 

Wolff, B., Koester, D., Dreizler, S., & Haas, S. 1998, A&A, 329, 1045 

Wolff, B., Koester, D., & Lallement, R. 1999, A&A, 346, 969 

Wolff, B., Kruk, J. W., Koester, D., AUard, N. F., Ferlet, R., & Vidal-Madjar, A. 2001, 
A&A, 373, 674 

Wood, B. E., Linsky, J. L., Hebrard, G., Vidal-Madjar, A., Lemoine, M., Moos, H. W., 
Sembach, K. R., & Jenkins, E. B. 2002, ApJS, 140, 91 

York, D. G. 1983, ApJ, 264, 172 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-27- 




910 



920 



930 



940 



950 




990 



LiF lA MORS 



1020 



1040 



1060 



1080 



12 









LiF 2A LWRS 



1100 



1180 



1120 1140 1160 

Wavelength (A) 

Fig. 1. — FUSE spectra of GD 246 with identification of interstellar lines. An annotation on 
each panel indicates which channels are displayed and if the data was obtained through the 
LWRS or MDRS apertures. Photosphcric lines are marked with dotted lines and geocoronal 
emission is annotated with 0. Ar I A1066 is blended with stellar Si IV. The data are binned 
by 4 for display purposes only. 
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Fig. 2. — Same as Fig. 1 but for WD 2331—475. The data presented here are from the 
MORS observation, P1044203. 
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Fig. 3.— Same as Fig. 1 but for HZ 21. 
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Fig. 4. — Same as Fig. 1 but for Lan 23. Many lines of molecular hydrogen (not marked) 
are also seen along this hne of sight, for example in the range 1000 - 1020 A. 
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Fig. 5. — STIS echelle spectra of GD 246, showing two absorption components. The two 
components at Vq = —9.8 and —0.6 km s^^ are marked by dashed hnes. The weak component 
at —0.6 km s~^ is identified with the LIC in §7.5. The bold hne represents the fit to the 
absorption hnes. 
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Fig. 6. — GD 246 - Fits to the D I lines used in the analysis. L and M stand for LWRS 
and MDRS data, respectively. D I Ly(3 was not used in the analysis due to contamination 
by geocoronal emission. Dashed vertical lines mark the centroids of the fitted D I and O I 
absorption lines. 
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Fig. 7. — Top panel: Interstellar H I Lya toward GD 246 (data binned by 4 for display 
purposes only). The histogram shows the coadded STIS E140M observations of GD 246. 
Model a is the best fit stellar model with Teg = 53,000 K and log g = 7.85 (sec text for 
discussion). Model b was computed with T^s — 50,000 K and log g = 7.7, and model c with 
Tcff = 56,000 K and log g = 8.0. Bottom panel: Fit to the H 1 Lya profile after division 
by the best fit stellar model (model a, see text) and with the continuum modeled by a 2"*^ 
order polynomial. 
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Fig. 8. — WD 2331—475 - Double component fits to some metal lines used in the analysis. 
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from the profile fitting and curve of growth analysis. In some cases, the discrepancy between 
the single component fit and the data was apparent upon inspection of the blue side of the 
absorption fine. Dashed vertical fines mark the centroids of the two components for C II*, 
N I, and Fe II. 
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Fig. 9.— WD 2331-475 - Fits to the D I lines used in the analysis, all from the P1044203 
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Fig. 11. — HZ 21 - Fits to the D I lines used in the analysis, a, b, and c represent data from 
P204 MDRS, P204 LWRS, and M108 LWRS, respectively. Dashed vertical hnes mark the 
centroids of the fitted D I and O I absorption lines. 
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Fig. 12. — Lan 23 - Fits to the D I lines used in the analysis. Dashed vertical hnes mark the 
centroids of the fitted D I, O I, and H2 absorption lines. 
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Table 1. Stellar properties 



Quantity 


GD 216 


AVD 2331-175 


HZ 21 


Lau 23 


WD number 


2309+105 


2331-475 


1211+332 


2247+583 


Spectral Type 


DA 


DA 


DO 


DA 


I (deg) 


87.25 


334.85 


175.04 


107.64 


b (deg) 


-45.11 


-64.81 


+80.03 


-0.64 


d" (pc) 


79 ± 24 


82 ± 25 


115 ± 35 


122 ± 37 


V 


13.09 


13.1 


14.22 


14.26 


Teff (K)^ 


53,088 ± 968 


51,800 ± 800 


53,000 ± 2650 


59,360 ± 800 


log g (cm s"^)'' 


7.85 ± 0.07 


7.79 ± 0.07 


7.8 ± 0.3 


7.84 ± 0.05 


Reference 


1 


2 


3 


2 



*A11 distances are photometric, 30% error adopted. 
''The uncertainties quoted are 1 a. 

References. — (1) Napiwotzki et al. (1999); (2) Venues et al. (1997); (3) Dreizler & 
Werner (1996) 
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Table 2. Log of FUSE observations 



Star 




Datp 






Flicn Timp 


Mode** 


S/N'^ 












fksoc) 






GD 246 


M1010602 


1999 Dec 09 


MDRS 


2 


1.2 


HIST 


4.6 




M1010604 


1999 Dec 10 


LWRS 


7 


3.4 


HIST 


6.4 




P1044101 


2000 Jul 19 


LWRS 


28 


14.8 


HIST 


8.6 




M1010601 


2000 Nov 12 


LWRS 


3 


1.6 


HIST 


5.2 




P2042401 


2001 Jul 14 


MDRS 


51 


24.6 


HIST 


8.3 


WD 2331-475 


P 1044202 


1999 Nov 08 


LWRS 


29 


20.2 


HIST 


8.2 




P1044201 


2000 Jun 23 


LWRS 


20 


19.4 


HIST 


8.4 




P1044203 


2002 Jul 27 


MDRS 


72 


32.1 


HIST 


8.3 


HZ 21 


M1080201 


2000 May 10 


LWRS 


5 


4.9 


TTAG 


4.4 




P2040801 


2001 Jan 27 


LWRS 


4 


12.4 


TTAG 


5.2 




P2040802 


2001 Feb 6 


MDRS 


25 


16.7 


TTAG 


5.5 


Lan 23 


P2510101 


2001 Jul 19 


LWRS 


8 


20.4 


TTAG 


4.7 



^See text for definition of aperture. 

^Scc text for definition of observing mode. 

'^S/N ratios, per detector pixel, measured in the LiF lA segment, between 1000-1005 A for all stars 
except Lan 23, where we use the range 995 - 1000 A. 
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Table 3. STIS echelle observations of GD 246 



Program ID 


Date 


Grating 


Aperture 


Exp. Time 


Wavelengths Covered 










(s) 


(A) 


O4G102010 


1998 Nov 20 


E140M 


0'.'02 X 0'/06 


1699.0 


1140-1735 


O4G102020 


1998 Nov 20 


E140H 


0'.'02 X 0'/09 


2420.0 


1170-1372 
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Table 4. Atomic data and analysis methods for the hnes used in the analyses" 



Species 


Wavelength (A) 


Log /A 


GD 246 


WD 2331-475 


HZ 21 


Lan 23 


H I 


1215.670 


2.70 


P 








D I 


920.713 


0.17 








P 




922.899 


0.31 








P 




925.974 


0.47 


P 






P 




930.495 


0.65 


P 




P 


P 




937.548 


0.86 


P 


P 


P 


P 




949.485 


1.12 


P 


P 


P 






972.272 


1.45 


P 


P 








1025.443 


1.91 




P 






C II* 


1037.018 


2.11 


P 


P, A 


A 


A 




1335.663 


1.23 


P 










1335.708 


2.23 


P 








N I 


952.303 


0.25 








P 




952.415 


0.21 








P 




952.523 


-0.24 








P 




953.415 


1.10 




P, C 


P, C 






953.655 


1.38 




P, C 


P, C 






953.9699 


1.52 




P, C 








954.1042 


0.81 




P 








963.990 


1.54 




P 


P, C 






964.626 


0.96 


P, C 


P, c 


P, C 






965.041 


0.59 


P, c 


p 


P, C 






1134.1653 


1.69 




P, c 








1134.4149 


1.53 




P, c 








1134.9803 


1.24 




P, c 








1199.550 


2.19 


c 










1200.223 


2.01 


c 










1200.710 


1.71 


c 








N II 


1083.994 


2.10 


A 


A 


A 


A 


I 


919.658 


-0.06 


P, C 


P, c 


P, C 






919.917 


-0.79 


P, c 


p 


P, C 


P, A 
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Table 4 — Continued 



Species 


Wavelength (A) 


Log /A 


GD 246 


WD 2331-475 


HZ 21 


Lan 23 




921.875 


0.04 


P, C 


P, C 








922.200 


-0.65 


P 


P 








924.950 


0.15 


P, C 


P, C 


P, C 






925.446 


-0.49 


P, c 


P 


P, c 






929.5168 


0.32 




P, c 








930.257 


-0.30 


P, c 


p 


P, c 






936.6295 


0.53 




P, c 








948.6855 


0.77 


c 


c 


c 


C 




950.885 


0.18 


P, c 










971.738 


1.13 


c 


c 


c 


C 




972.142 


-0.47 


p 


p 








974.070 


-1.82 








P 




976.4481 


0.51 




P, c 








1039.2301 


0.98 


c 


c 


c 


C 


Si II 


1020.699 


1.22 


P, c 


P, A 


P, A 


A 




1190.416 


2.54 


P, c 










1193.230 


2.84 


P, c 










1 on /I QTn 
ioU4.o lu 


o no 
z.Uo 


P, L 








P II 


963.800 


3.15 


P 


P, A 


P, A 


A 


S II 


1250.584 


0.83 


P, c 










1253.811 


1.14 


P, c 










1259.519 


1.32 


P, c 








Ar I 


1048.220 


2.44 


p 


P, A 


P, A 


A 


Fell 


1063.176 


1.76 


p 




P, A 


P, C 




1121.975 


1.36 








P, C 




1125.448 


1.26 








P, C,A 




1144.938 


2.08 


p 


P, A 




P, C 



^P, C, and A, denote lines that are analyzed with profile fitting, curve of growth, and 
apparent optical depth, respectively. 
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Table 5. Adopted column densities' 



Species 


GD 246 (cm-2) 


WD 2331- 


-475 (cm"^) 


HZ 21 


(cm ^) 


Lan 23* (cm'^) 


log N(H I) 


19.11 


± 


0.05 












19.89 ± 


log N(D I) 


14.29 


± 


0.09 


14.19 


± 0.12 


14.40 


± 


0.15 


15.23 ± 0.13 


log N(C II*) 


13.05 


± 


0.04 


13.24 


± 0.17 


12.97 


± 


0.20 
0.27 


> 13.65 


log N(N I) 


14.75 


± 


0.06 


14.53 


± 0.10 


14.77 


± 


0.08 


15.73 ± ^il 


log N(N 11) 


> 


13.J 


Bl 


> 


14.27 


> 


14.03 


> 13.95 


log N(0 1) 


15.67 


± 


0.07 


15.48 


± 0.11 


15.74 


± 


0.10 


16.72 ± °:i 


log N(Si 11) 


14.07 


± 


0.06 
0.04 


13.99 


± 0.20 


14.36 


± 


0.08 


> 14.52 


log N(P 11) 


12.29 


± 


0.10 


12.18 


± 0.20 


12.57 


± 


0.12 
0.15 


> 12.72 


log N(S 11) 


14.34 


± 


0.02 














log N(Ar I) 


13.14 


± 


0.13 
0.10 


12.77 


± 0.20 


13.16 


± 


0.05 


> 13.50 


log N(Fe II) 


13.30 


± 


0.10 


13.32 


± 0.10 


13.54 


± 


0.08 


14.03 ± 0.13 



^Adopted values are the combination of PFl and PF2 (see text for more on this). All 
uncertainties are 2a. 

'^N(H I) from Wolff et al. (1999). Range quoted includes uncertainty in the photospheric 
composition of this star (see §5). 
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Table 6. H2 column densities in Lan 23" 



logN(H2) Profile Fitting'^ Curve of Growth^ 



J = 





14.35 ± 0:°^ 


14 12 ± 

-'-^•-'-^ 0.09 


J = 


1 


14.97 ± ',^1 


14.89 ± 0.11 


J = 


2 


14.32 ± HI 


14.20 ± ^1 


J = 


3 


13.95 ± 0.07 


13.93 ± q;^q 



^AU uncertainties are 2a. 

^The b value of the fit was constrained to 
4.1 km s~^. The errors quoted do not include a 
contribution from errors in the b value. 

^All hnes are fit simultaneously yielding 
6 = 4.lt°| km s-i (2a). 
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Table 7. Ratios of column densities" 



Quantity GD 246 WD 2331-475 HZ 21 Lan 23 



D I/H I (xlO- 


-5-\ 


1.51 


± 


0.39 
0.33 














2.19 


± 


4.64 
0.68 


I/H I (xlO- 


-4\ 


3.63 


± 


0.77 
0.67 














6.76 


± 


15.5 
4.10 


D I/O I (xlO- 


-2\ 


4.17 


± 


1.21 
1.00 


5.13 


± 


2.20 
1.69 


4.57 


± 


2.22 
1.63 


3.24 


± 


3.27 
2.06 


N I/H I (xlO- 


-5\ 


4.37 


± 


0.84 
0.74 














6.92 


± 


14.7 

2.33 


D I/N I (xlO- 


-1\ 


3.47 


± 


0.95 
0.79 


4.57 


± 


1.88 
1.45 


4.27 


± 


1.96 
1.44 


3.16 


± 


1.56 
1.23 


I/N I 




8.32 


± 


1.91 
1.64 


8.91 


± 


3.45 
2.71 


9.33 


± 


3.07 
2.48 


9.77 


± 


9.89 
6.37 



^AU uncertainties are 2a. 
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Table 8. Revised FUSE ratios and comparison with other measurements" 



Quantity Moos ct al. (2002) for Degrees of Meyer ct al. 

+ this work'' mean freedom v'^ (1997, 1998)"^ 



D I/H I 


1.52 


± 


0.07 (0%) 


1.0 


5(4) 








(xlO-5) 


















I/H I 


3.17 


± 


0.19 (+5%) 


3.5 


5(4) 


3.43 


± 


0.15 


(xlO-4) 


















D I/O I 


4.06 


± 


0.17 (+2%) 


1.9 


9(6) 








(xlO-2) 


















N I/H I 


4.28 


± 


0.25 (+1%) 


1.4 


4(3) 


7.5 


± 


0.4 


(xlO-5) 


















D I/N I 


3.41 


± 


0.15 (+1%) 


2.2 


9(5) 








(xlO-i) 


















I/N I 


8.23 


± 


0.38 (+2%) 


2.9 


8(5) 


4.6 


± 


0.3 



^All uncertainties in this table are la in the mean. 

^Mean values and uncertainties in the mean. The values in () corre- 
spond to the percent increase from the Moos et al. (2002) values. Our 
values are combined with those from Moos et al. (2002) by taking a 
weighted mean. We use the largest of the lower and upper error bars of 
each individual ratio to compute the weighted mean. 

'^y — number of sightlines — 1. Values in () represent the number of 
degrees of freedom in Moos et al. (2002) 

•^The values in the table correspond to 0/H and N/H, the uncertainties 
quoted are errors in the mean. 0/N calculated here from 0/H and N/H. 
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Table 9. Comparison of abTindanccs toward GD 246 with solar abundances 



X log(X/H)Q +12 log(X/H)^^ 246 Reference 



N 


7.93 


7.64 ± I f 


1 





8.69 


8.56 ± HI 


2 


Si 


7.54 


6.96 ± \il 


1 


p 


5.51'' 


5-18 ± \il 


3 


s 


7.27^ 


7.23 ± \Z 


3 


Ar 


6.40 


6.03 ± \il 


3 


Fe 


7.45 


6.19 ± lit 


1 



^All uncertainties are 2cr. 

^The meteoritic and photospheric values were aver- 
aged. 



References. — (1) Holweger (2001); (2) AUende Pri- 
eto et al. (2001); (3) Grevesse & Sauval (1998). 



